Reactions of ultracold alkali metal dimers 
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We investigate the energetics of reactions involving pairs of alkali metal dimers. Atom exchange 
reactions to form homonuclear dimers are energetically allowed for some but not all of the het- 
eronuclear dimers. We carry out high-level electronic structure calculations on the potential energy 
surfaces of all the heteronuclear alkali metal trimers and show that trimer formation reactions are 
always energetically forbidden for low-lying singlet states of the dimers. The results have important 
implications for the stability of quantum gases of alkali metal dimers. 

PACS numbers: 34.20.-b, 34.50.Cx, 37.10.Pq 



It has recently become possible to create samples of al- 
kali metal dimers in deeply bound states at temperatures 
below 10-6 K [H-i]. For KRb \^ and Csa 1^], molecules 
are first formed in high-lying vibrational states by magne- 
toassociation and then coherently transferred to the abso- 
lute ground state by stimulated Raman adiabatic passage 
(STIRAP). These capabilities open up new possibilities 
for ultracold chemistry, for creating strongly interacting 
quantum gases, and for producing tunable models of im- 
portant condensed-phase systems [loj . 

Ultracold molecules are usually confined in very shal- 
low traps. Any collision that transfers internal energy 
into relative kinetic energy is likely to eject both colli- 
sion partners from the trap. If both species are in their 
absolute ground state, inelastic collisions are impossible, 
but there remains the possibility of reactive collisions. 
Indeed, for fermionic *°K^^Rb, Ospelkaus et al. [lH have 
carried out detailed studies of the exothermic atom ex- 
change reaction, 



KRb + KRb 



Ko + Rb,. 



(1) 



When all the molecules are in the same nuclear spin state, 
the reaction rate is strongly suppressed by the Pauli prin- 
ciple. However, if some of the molecules are transferred 
into a different spin state, the reaction proceeds very fast 
and the molecules are lost from the trap. 

As will be seen below, atom exchange reactions analo- 
gous to (HI) are energetically allowed for some alkali metal 
dimers but forbidden for others. However, even when 
atom exchange is forbidden, there remains in principle 
the possibility of atom transfer reactions such as 



KRb + KRb — ^ K + KRba 



KaRb + Rb. (2) 



In a simple pairwise-additive model of the energet- 
ics, the reactants in such a reaction have two nearest- 
neighbor interactions and the products have three, so 
some researchers have anticipated that the trimer forma- 
tion reactions would be energetically allowed. However, 
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pairwise-additive models are known to be very poo r for 
the quartet excited states of alkali metal trimers [IJ] and 
are likely to be even poorer for the doublet ground states. 
The principle purpose of this paper is to explore the ener- 
getics of trimer formation reactions such as ^ . We will 
demonstrate that, for singlet alkali metal dimers in levels 
near the potential minimum, trimer formation reactions 
are in fact always energetically forbidden. 

Before proceeding to the trimer formation reactions, 
we briefly consider atom exchange reactions analogous to 
^ for the heteronuclear dimers formed from the alkali 
metals Li, Na, K, Rb and Cs. All the homonuclear and 
heteronuclear dimers except LiRb have been studied in 
detail by high-resolution spectroscopy, and dissociation 
energies accurate to ±1 cm"^ or better have been 
extracted as listed in Table HI The energy changes for 
the atom exchange reactions can therefore be calculated 
directly from experiment, and are summarized in Table 
im The values given are taken from dissociation energies 
Dg measured to the dimer equilibrium geometries and 
so are subject to small corrections for the differences in 
zero-point energy between reactants and products. These 
corrections can be up to -f 25 cm"^ for LiX systems but 
are less than ±2 cm^^ for the remainder. It may there- 
fore be concluded that all the heteronuclear Li dimers and 
KRb will be subject to reactive trap loss, but all the re- 
mainder should be stable with respect to atom exchange 
collisions in their ground rovibronic state. 

Trimer formation reactions cannot be considered in a 
similar way because an experimental binding energy is 
available only for Lis [23] and not for any of the het- 
eronuclear trimers. We have therefore carried out elec- 
tronic structure calculations for all the homonuclear and 
heteronuclear alkali metal trimers, using the multirefer- 
ence average-quadratic coupled-cluster method (AQCC). 
All calculations used the MOLPRO package [2^. The 
alkali atoms were described in a single-electron model 
and the core- valence interaction was taken into account 
using an effective core potential (ECP) with a core po- 
larization potential (CPP). We used the ECPxSDF fam- 
ily of core potentials, developed by the Stuttgart group 
[29I [30I I , with core polarization potentials based on those 
of Miiller and Meyer [3l| . We obtained modified values 
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TABLE I: Dissociation energies De (in cm ^) for alkali metal dimers. The quantities in parentheses are uncertainties in the 
final digit(s). 



Li 

Na 

K 

Rb 

Cs 



Li 

8516.768(8)" 



Na 

7105.5(1.0)' 
6022.0286(53)^ 



K 

6216.886(100)" 
5273.62(10)" 
4450.906(50)^ 



Rb 
5946(100)"^ 
5030.502(10)'' 
4217.815(10)* 
3993.47(18)™ 



''Ref.lil 
'^Ref. [l| 

"^The binding energy for LiRb is not available from experiment, so 
this value is calculated using the AQCC method described in this 
paper. 
'Ref.Iii 
/Ref. g 



SRef. 
''Ref. 
'Ref. 
JRef. 
^Ref.m 
'Ref. 
"Ref. 
"Ref.L. 

°Do from Ref. and zero-point energy from Ref. [26 



Cs 

5875.542(5)' 
4954.237(100)' 
4069.208(40)' 
3836.14(50)" 
3649.695(2)° 



TABLE II: Energy changes /S.E2 for the reactions 
2XY— >X2+Y2 (in cm~^). The quantities in parentheses are 
uncertainties in the final digit (s). 





Na 


K 


Rb 


Cs 


Li 


-328(2) 


-533.9(3) 


-618(200) 


-415.38(2) 


Na 




74.3(3) 


45.5(5) 


236.75(20) 


K 






-8.7(9) 


37.81(13) 


Rb 








29.1(1.5) 



of the Miiller-Meyer cutoff parameter (0.95 for Li, 0.82 
for Na, 0.36 for K, 0.265 for Rb and 0.24 for Cs) that 
reproduce the experimental bond lengths of the ground- 
state homonuclear alkali dimers at the same level of the- 
ory. We used the uncontracted s» basis sets designed for 
ECPxSDF core potentials [H, IsOl. augmented by addi- 
tional s, p, d and / functions [33 • With these polar- 
ization potentials and basis sets we reproduced the sin- 
glet binding energies for homonuclear alkali metal dimers 
with an accuracy better than 1% for Li2, Na2, Rb2 and 
2% for K2 and Cs2. The binding energies for the het- 
eronuclear dimers are as good as for the homonuclear 
dimers, except for LiCs, for which the error in the bind- 
ing energy was -1-2.5%. 

To understand the doublet states of heteronuclear 
alkali metal trimers, it is useful first to consider the 
homonuclear systems. The important molecular orbitals 
are those formed from the outermost s orbitals on each 
atom. At an equilateral triangular configuration {D^h 
symmetry), the two lowest-lying molecular orbitals of a 



TABLE III: Atomization energies and equilibrium geometries 
of the X2Y trimers from AQCC calculations, together with 
energy changes Ai?3 for the reactions 2XY— >-X2Y + Y, ob- 
tained by combining the trimer results with dimer dissociation 
energies obtained with the same method. 



Atomization 


energy ( 


cm-^) 




X 






Li 


Na 


K 


Rb 


Cs 


Li 


13189 


9977 


8341 


7982 


8378 


Na 


11583 


8113 


7125 


6783 


7140 


Y K 


10681 


7795 


6258 


5902 


5890 


Rb 


10499 


7649 


6080 


5685 


5661 


Cs 


11073 


8128 


6211 


5781 


5494 



r-XiY.rxaY.rxx (A) X 

Li Na K Rb Cs 

Li 2.8,2.8,3.2 3.0,3.0,4.0 3.5,3.5,5.3 3.6,3.6,5.9 4.1,4.1,4.7 

Na 3.1,3.5,2.7 3.2,3.2,4.2 3.7,4.4,4.0 4.0,4.4,4.2 4.1,4.4,4.6 

Y K 3.5,4.3,2.8 3.7,3.7,3.9 4.1,4.1,5.2 4.2,4.2,5.7 4.4,5.5,4.8 

Rb 3.6,4.5,2.8 3.8,3.8,3.8 4.2,5.3,4.1 4.4,4.4,5.5 4.6,5.5,4.7 

Cs 3.8,3.8,3.1 4.0,4.0,3.7 4.5,4.5,4.9 4.6,4.6,5.5 4.8,4.8,5.7 



AEz (cm-i) 






X 








Li 


Na 


K 


Rb 


Cs 


Li 


3759 


4145 


3979 


3910 


3660 


Na 


2539 


3843 


3281 


3287 


2962 


Y K 


1639 


2611 


2460 


2444 


2264 


Rb 


1393 


2421 


2266 


2295 


2101 


Cs 


965 


1974 


1943 


1981 


1958 
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FIG. 1: (Color online) Upper panel: The intersecting poten- 
tial energy surfaces for A\ and states of Rb2Cs at C^v 
geometries, as a function of the RbCs bond length r and the 
bond angle Q. Lower panel: curves for the A\ and ^B^ states 
of Rb2X systems as a function of bond angle Q, with the bond 
length optimized at each angle. For X=Li,Cs,K the minima 
for ^B^ states are the absolute minima, while for Rb2Na sys- 
tem the absolute minimum originates from the distorted A\ 
minimum geometry. 



homonuclear trimer have a\ and e symmetry. The low- 
est doublet state has configuration a{e^ . It is therefore 
orbitally degenerate, with '^E symmetry, and is subject 
to a Jahn- Teller distortion to an isosceles geometry {Civ) 
that splits the e orbitals into a\ and 62 components: the 
62 orbital has a node between the two equivalent atoms. 
The equilibrium structures of the homonuclear trimers 
all have Civ geometries with ground states of sym- 
metry. 

For a heteronuclear trimer X2Y, the symmetry is al- 
ways Civ or Cg. For Civ geometries, the upper a\ orbital 
and the 62 orbital are close together in energy and the 
minimum (restricted to Civ) may be on either the A\ 
surface or the ^Bi surface. We have therefore calculated 
the energy for all the heteronuclear trimers in both A\ 



and ^Bi states for Civ geometries. Typical results are 
shown for Rb2Cs in the upper panel of Fig. [TJ The ge- 
ometry is specified by a bond length r = rxjY ~ ^X2Y 
and the angle d between the two XY bonds. It may 
be seen that the two surfaces intersect at an angle near 
Q — 50°: since the two states have the same symmetry at 
Cs geometries, they actually intersect only at Civ geome- 
tries, producing a seam of conical intersections there. An 
alternative representation of the results, for all the Rb2X 
systems, is shown in the lower panel of the Figure: in 
this case r has been optimized to find the energy mini- 
mum for each value of 0, producing intersecting potential 
curves rather than 2-dimensional surfaces. The minima 
on the two surfaces are usually close together in energy 
(always within 1000 cm~^, but often within 200 cm~^). 
However, the ^Bi minimum is below the A\ minimum 
for all the trimers except the 7 heteronuclear X2Na and 
CS2X species; for Rb2Cs, shown in Fig. [1] the ^Bi min- 
imum is near Q = 63°. The equilibrium geometries and 
energies for both states are provided as supplementary 
material [33| . 

For heteronuclear trimers there is the additional pos- 
sibility of distortion to a lower-symmetry Cs (scalene) 
geometry. We have therefore explored whether such dis- 
tortions lower the trimer energies. At Cs geometries the 
valence orbitals formed from atomic s orbitals are all of 
a! symmetry, so both low-lying states have A! symme- 
try and can mix. Nevertheless, in most cases it is clear 
whether the singly-occupied orbital has bonding charac- 
ter (fli-like) or antibonding character (62-hke) between 
the two like atoms. For Cs2Li, where the A\ state was 
already below the ^Bi state, distortion does not lower the 
energy and the equilibrium geometry has Civ symmetry. 
However, for all the other systems the geometry corre- 
sponding to the A\ minimum is in fact a saddle point on 
the full 3-dimensional surface: for Li2Na, K2Na, Rb2Na, 
Cs2Li, Cs2Na, CS2K and Cs2Rb, this simply deepens the 
minimum. For Li2K, Li2Rb, K2Rb, K2CS, Rb2Na, the 
distortion produces a A! state whose absolute minimum 
(of Cs symmetry) is lower in energy than the ^Bi state 
(which always retains an equilibrium geometry of Civ 
symmetry). However, for Li2Cs, Na2Li, Na2K, Na2Rb, 
Na2Cs, K2Li, Rb2Li, Rb2K and Rb2Cs the energy gained 
by distortion is not enough and the "^Bi state of Civ ge- 
ometry remains the absolute minimum. 

Table [TTTl summarizes the trimer atomization energies, 
equilibrium geometries and the energy change for the 
trimer formation reactions for all the alkali metal trimers 
from Li to Cs. It may be seen that all the trimer for- 
mation reactions (from singlet dimers) are substantially 
endoergic. Trimer formation reactions will therefore not 
take place for alkali metal dimers formed in singlet states 
near the bottom of the potential well. However, trimers 
may of course still be formed from dimers in triplet states, 
which are much more weakly bound, or from high-lying 
vibrational states, including Feshbach molecules. 

The trimer energies are always substantially below the 
energy of any atom -t- diatom arrangement of the same 
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atoms. The entrance channels of chemical reactions be- 
tween alkali metal atoms and singlet dimers are thus 
likely to be barrierless, as shown by Tscherbul et at. [s^l 
for Rb + RbCs (though of course the reactions them- 
selves will not always be energetically allowed). However, 
a full treatment of the dynamics of these reactions would 
require a detailed study of the potential energy surfaces 
for at least the two lowest-lying electronic states and the 
interactions between them. This contrasts with the situa- 
tion for reactions involving spin-stretched states of alkali 
metal atoms and triplet dimers, which have been studied 
using single electronic surfaces for the quartet states of 
the trimers [H-lil]. 

Our atomization energies for homonuclear systems 
may be compared with 13436 cm^^ for Lia from multi- 
reference configuration interaction (MRCI) calculations 
[jlj . and 5437.1 cm^^ for Csa from full configuration in- 
teraction (CI) calculations [42]. Our values for trimers 
containing Li, Na and K also agree well (within 1000 
cm~^) with early CI work by Pavolini and Spiegelmann 
[43l |. In all cases the calculations used effective core po- 
tentials similar to those in the present work. 

The present results for trimer energies may be ratio- 
nalized using a very simple model. In the simplest form 



of Hiickel theory, considering only one s orbital on each 
atom, with a bond integral /3, a homonuclear dimer in 
a singlet state has binding energy 2|/3|. An equilateral 
trimer has binding energy 3|/3|, while a linear trimer has 
binding energy 2-\/2|/3|. An atom transfer reaction such 
as is therefore endoergic by \f3\ or slightly more, i.e. 
by about half the dimer binding energy. This is quite 
different from the result predicted by pairwise additivity. 
However, simple orbital-based models of chemical bond- 
ing must be treated with caution for the alkali metals, 
because they have low-lying p orbitals that often make 
important contributions to bonding. Ion-pair states can 
also be important. Thus, while Hiickel theory can be used 
to rationalize the results of the present work, it could not 
have been used to predict them. The use of high-level 
electronic structure calculations, as in the present work, 
is essential to obtain reliable conclusions. 
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